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ABSTRACT 

Accretion onto compact objects plays a central role in high energy astrophysics. In 
these environments, both general relativistic and plasma effects may have significant 
impacts upon the spectral and polarimetric properties of the accretion flow. In paper 
I we presented a fully general relativistic magnetoionic theory, capable of tracing rays 
in the geometric optics approximation through a magnetised plasma in the vicinity of 
a compact object. In this paper we discuss how to perform polarised radiative transfer 
along these rays. In addition we apply the formalism to a barotropic thick disk model, 
appropriate for low luminosity active galactic nuclei. We find that it is possible to 
generate large fractional polarisations over the innermost portions of the accretion 
flow, even when the emission mechanism is unpolarised. This has implications for 
accreting systems ranging from pulsars and X-ray binaries to AGN. 
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1 INTRODUCTION 

The spectral and polarimetric properties of astrophysical ob- 
jects can provide significant insights into their structure and 
dynamics. As a result, a number of theoretical investigations 
into the source of these properties have been undertaken. 
Many of these have been primarily concerned with the spec- 
tral properties alone, typically comparing a physically mo- 
tivated accretion flow to observations. However, with the 
measurement of polarisation in a number of sources, a sig- 
nificant fraction of the focus has been turned towards re- 
producing their polarimetric properties. In the context of 
an accreting compact object, both general relativistic and 
plasma effects can play a role in det ermining these proper- 
ties. In lBroderick fc Blandfor'dll2003l (hereafter Paper I) we 
demonstrated how to construct ray trajectories, in the geo- 
metric optics approximation, in a magnetoactive plasma in 
a relativistic environment. In order to apply this to realistic 
accretion environments it is necessary to be able to perform 
radiative transfer along these rays. 

Non-refractive, polarised radiative transfer through 
magnetised plasmas is flat space has been extensively stud- 
ied. A number of examples involving weak magnetic fields 
exist in the lit e rature (see e.g. IS azonov & Tsvtov ichlll968t 
ISazonovl Il969l: iJones fc 0'DelllTl977allbl: IGinzburJfigTI) . 
More recently, investigations into the net effects of tangled 
magnetic fields (expected to be t ypical in magnetised ac- 
cretion flows) have begun (see e.g. lRuszkowski fc Beeelmanl 
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120021) . However, none of these deal with general relativistic 
environments. 

The importance of refraction in the propagation of ra- 
dio wavelengths has lo ng been appreciated in the con text of 
the ionosphere (see e.g. lBudderll964l : lGinzbur Jl97Ci) . More 
recently, refraction has been stud ied in conjunction with pul- 
sars (see e.g. [Wcltevrede et aljEp 03; Pctrova 2002, 20o3 
iBarnard fc Arondlll986l: lArons fc Barnardlll98Q') . Nonethe- 
less, in all of these cases, the emission was assumed to orig- 
inate from a region distinct from where the refraction oc- 
curred. Re fractive lensing of neutron stars was considered by 
ih aviv et al. 1999, but ignored general relativisitic effects. 

General relativistic studies into the propagation of po- 
larisation in vacuum have been done. These have been 
primarily interested in the geometrical effects due t o the 
parall el transport of the linear p olarisation (see e.g. lAgoll 
.1997tlLaor et al.lll990l; IConnors et al. .1980ii) . Alternatively, 
iir iBromle^^^nl20ML polarised emission in a general rel- 
ativistic environment is considered. However, none of the 
typical plasma transfer effects (e. g. Faraday rotat ion) were 
included in these calculations. In lHevl et al . 2003, the vac- 
uum birefringence due to strong magnetic fields was con- 
sidered in the context of neutron star atmospheres. How- 
ever, in both, refraction was completely ignored. There have 
been some attempts to study t he pro blem of ray propagation 
in a covariant form (see e.g. iBroderick fc BlandfordI 1200^ : 
iMelrose fc GedalinllioOlh . but in these the radiative trans- 
fer was not addressed. 

As discussed in Paper I, refraction coupled with the 
presence of a horizon can be a source of significant polar- 
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isation when the observation frequency is near the plasma 
and cyclotron frequencies of emitting region. The sense of 
the resulting net polarisation is determined by the plasma 
parameters at the surface at which the polarisation freezes 
out (when the modes cease to be adiabatic and must be 
treated as if they were in vacuum). Typically, this will re- 
sult in a net circular polarisation. In a future paper we will 
discuss astrophysical environments in which this may be the 
case, including applications to Sgr A*and high mass X-ray 
binaries. 

We present a method for performing polarised radiative 
transfer through a strongly refractive magnetised plasma in 
a general relativistic environment. Additionally, we apply 
this to a model of a thick accretion disk. This is done in six 
sections with J^lbriefly reviewing the formalism presented in 
Paper I, f|3discussing how to perform the radiative transfer 
in a magnetised plasma, ^1] presenting low harmonic syn- 
chrotron as a possible emission mechanism, JS] presenting 
some results, and containing conclusions. The details of 
constructing a magnetised, thick, barotropic disk are pre- 
sented in the appendix. 



2 RAY PROPAGATION 

While astrophysical plasmas will, in general, be hot, the cold 
case provides an instructive setting in which to demonstrate 
the types of effects that may be present. As a result, it will 
be assumed that the plasma through which the rays prop- 
agate will be cold, with a small component of emitting hot 
electrons. As shown in Paper I, the rays may be explicitly 
constructed given a dispersion relation, D{k^,x'^) (a func- 
tion of the wave four-vector and position which vanishes 
along the ray), by integrating the ray equations: 



"d7 



dP 
dkn 



and 



"d7 



dP 



(1) 



where r is an afhne parameter along the ray. Expand- 
ing Maxwell's equations in the geometric optics limit pro- 
vides the polarisation eigenmodes and the dispersion rela- 
tion (given a conductivity): 



{k^k^SH: - k^k^ - ATviioa"^) E" = , 



(2) 



where i?^ is the four-vector coincident with the electric 
field in the locally flat, comoving rest frame (LFCR frame), 
ui = -~u'^kf_i {u'^ is the plasma four- velocity which defines the 
LFCR frame), and cr^^ is the covariant extension of the the 
conductivity tensor. For the cold, magnetoactive, electron- 
ion plasma (in the limit of infinite ion mass), the dispersion 
relation is 



Pik^^x") 



k'^kn, — Suj 



2(1 + 5) 



\ mijj 



- (1 + 25) ujI 



\ mui 



(3) 



where B'^ is the four-vector coincident with external mag- 
netic field in the LFCR frame, a;p is the plasma frequency 



in the LFCR frame, ujb is the cyclotron frequency associ- 
ated with B'^, and 5 = uj%/{uj% — cop), This is a covariant 
form of the Applet on-Ifartree dispersion relation (see e.g. 
iBovd fc Sandersonl flQeQl . 

In general, the electromagnetic polarisation eigenmodes 
will not follow the same trajectories, and in particular will 
not follow null geodesies. As a result, the different polarisa- 
tion eigenmodes will sample different portions of the accre- 
tion flow. As shown in Paper I, it is possible for one mode 
to be captured by the central black hole while the other 
escapes, leading to a net polarisation. 



3 POLARISED RADIATIVE TRANSFER IN 
REFRACTIVE PLASMAS 

Both emission and absorption are local processes. However, 
because the transfer of radiation necessarily involves a com- 
parison between the state of the radiation at different points 
in space, global propagation effects need to be accounted 
for. These take two general forms: correcting for the gravi- 
tational redshift; and keeping track of the local coordinate 
system, i.e. ensuring that polarised emission is being added 
appropriately in the presence of a rotation of the coordinate 
system propagated along the ray. In addition, for a magne- 
toactive plasma, it is necessary to determine how to perform 
the radiative transfer in the presence of refraction. 



3.1 Length Scales and Regimes 

The problem of performing radiative transfer in a magne- 
toactive plasma has been treated in detail in the context of 
radio- wave pro pagation in the ionosphere (fo r a detailed dis- 
cussion see e.g. lGinzburJl97ol : lBuddenll964H . In these cases 
it was found that there were two distinct limiting regimes. 
These can be distinguished by comparing two fundamental 
scales of the afRne parameter r: that over which the polari- 
sation eigenmodes change appreciably, ts, and the Faraday 
rotation length, tf. Before rg can be defined it is necessary 
to define a pair of basis four-vectors that define the axes of 
the ellipse: 

lk°'k„+Lj'^)B'' -B^k^ik" -LJU") 

e(f = / (4) 



^kl^kfi + Lj2^(/fc-/fc<, + Lu^) B-iB^ - {B-ik^f 



e^""^l'u,k^Bn 



^ {k^ka + Lo'^) B-tB^ - {B-fk^f 



(5) 



where e'^'^"'' is the Levi-Civita pseudo-tensor. In terms of 
these, the ellipticity angle x can be defined by 



tanx 



^ eY.Ex. 



(6) 



In general, an additional angle, 0, is necessary to define the 
polarisation, namely the angle which defines the orientation 
of the ellipse. The basis four-vectors have been chosen such 
that (f) is identically zero. However, this choice introduces a 
new geometric term into the equations which accounts for 
the necessary rotation of the basis four-vectors, contribut- 
ing a non-zero dcjy/dr (see H3.3I for more details). Then, in 
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dx 
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J? 


dr 



'1/2 



(7) 



For the ordered fields employed here (see the appendices), 

(8) 



TS 



Lo'-^ dxi^ dr 



where this approximation form is true for small cyclotron 
and plasma frequencies and all but the most oblique angles 
of incidence. The Faraday rotation length is defined to be 
the distance over which the phase difference between the two 
polarisation eigenmodes reaches 2n, i.e. 



TF 



Afc, 



d^ 

' dT 



(9) 



where Ak^ is the difference between the wave vectors of the 
two modes. Strictly speaking in addition to tf, ts should 
be compared to a term describing the rate of change of the 
Faraday rotation length, however in the situations under 
consideration here this term is completely dominated by tf ■ 

Together, these length scales define three regimes: 
the adiabatic regime {tf <C Tg), the intermediate regime 
{2tf ~ Ts), and the strongly coupled regime {tf ^ ts). 
In all regimes the polarisation of the plasma eigenmodes is 
uniquely set by the dispersion equation, equation JJJ. 

In general, as 6 ^ 7r/2, AA; ~ (oipajs/cj^c) cos S + 
{uj%uj% /io^ c) , where 6 is the angle between the wave-vector 
and the magnetic field. Hence to remain in the adiabatic 
regime ts ^ {uj/ujb)'^tf{0 = 0), which is typically not true 
in astrophysical sources. As a result, as the magnetic field 
becomes perpendicular to the wave- vector, the modes gener- 
ally become strongly coupled. This is the reason why, when 
dealing with a large number of field reversals (e.g. in a molec- 
ular cloud) , the amount of Faraday rotation and conversion 
is oc B ■ dx and not jBj ■ dx (which would follow in the adi- 
abatic regime) despite the fact that Ts ^ tf{0 — 0) may be 
true throughout the entire region. 



3.2 Adiabatic Regime 

In the adiabatic regime the two polarisatio n modes propa- 
gate independently (see e.g. lGinzburdll97ol) . As a result, to 
a good approximation, the polarisation is simply given by 
the sum of the two polarisations. The intensities, lo and Ix, 
of the ordinary and the extraordinary modes, respectively, 
are not conserved along the ray due to the gravitational 
redshift. Consequently, the photon occupation numbers of 
the two modes. No and Nx, which are Lorentz scalars, and 
hence are conserved along the rays, are used. Therefore, the 
equation of radiative transfer is given by 

dA^o,x dl 



dr 



- d7 ^^O'^ 



ao,xNo,x) 



where 



d7 



dx'^ da;" 
dr dr 



dx>^ 



(10) 



(11) 



is the conversion from the line element in the LFCR frame to 
the affine parameterisation, and Jq x is the emissivity in the 
LFCR frame scaled appropriately for the occupation number 



(as opposed to the intensity). In practice, the occupation 
numbers will be large. However, up to fundamental physical 
constants, it is permissible to use a scaled version of the 
occupation numbers such that No,x = uj~^Io,x in vacuum. 

It is also this regime in which Faraday rotation and con- 
version occur. However, because these propagation effects re- 
sult directly from interference between the two modes, and 
hence require the emission to be coherent among the two 
modes, when they diverge sufficiently the modes must be 
added incoherently and thus Faraday rotation and conver- 
sion effectively cease. The modes will have divereged suffi- 
ciently when 



\Ax, 



AA 



(12) 



where AA is the emission band-width. For continum emis- 
sion, this reduces to lAx^l > A. Therefore in a highly refrac- 
tive medium an additional constraint is placed upon Faraday 
rotation. The depth at which equation (I12II is first satisfied 
can be estimated by considering an oblique ray entering a 
plane-parallel density and magnetic field distribution (at an- 
gle ( to the gradient). In this case, to linear order in ujp and 



As a result, 



\Axa 



d^Ax± 



UJBl-iJpZ 

2a;3 ' 



-smC— ~ — - — 
OZ uj-^z 



hence 



2aj3 
^ T 

UJBl^p 



(13) 



(14) 



The resulting number of Faraday rotations, wf, is then given 
by, 



np 



Ak 



-dz 



2-K 2-K sin 



(15) 



which is typically small for all but the smallest C^. Because, 
as discussed in section 13 linear polarisation is strongly sup- 
pressed by refraction, such a small Faraday rotation in neg- 
ligible. As a result, for the situations of interest here, in 
this regime the modes can be added together incoherently 
to yield the net polarisation. 



3.3 Strongly Coupled Regime 

In the limit of vanishing plasma density it is clear that the 
polarisation propagation must approach that in vacuum re- 
gardless of the magnetic field geometry. In this limit the two 
modes must be strongly coupled such that their sum evolves 
as in vacuum. In particular, it is necessary to keep track of 
their relative phases. This can be most easily accomplished 
by using the Stokes parameters to describe the radiation. In 
this case also it is possible to account for the gravitational 
redshift by using the photon occupation number instead of 
intensities. A'', Nq, Nu, Ny- However, it is also necessary to 
define the Nq, Nu, and A'^v' in a manner that is consistent 
along the entire ray. In order to do this we may align the 
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axes of Nq along the magnetic field, i.e. 

NQ=N{e';.)^N{el) 



Nu = N 



/ 1 1 - 
— =p.„ =f 



(16) 



where N{e'^) is the occupation number of photons in the 
polarisation defined by e'^ . Thus the problem of relating Nq , 
Nu, and Nv along the ray is reduced to propagating ej^ and 
e'^ . A change in r by dr is associated with a rotation of the 
basis by an angle 

dx" 



dr 



■Vue't.dT, 



(17) 



where the use of the covariant derivative, Vi,, accounts for 
the general relativistic rotations of ej^ and e^. As a result, 
the transfer effect due to general relativity and the rotation 
of the magnetic field about the propagation path is 



dNQ 
dr 

dNu 
dr 



(18) 



where the factor of 2 arises from the quadratic nature of N. 

After a specific emission model is chosen the emissivi- 
ties and the absorption coefficients are scaled as in H3.2I An 
example will be discussed in more detail in [21 

3.4 Intermediate Regime 

At some point it is necessary to transition from one limiting 
regime to the other. In this intermediate regime the polari- 
sation freezes out. A great deal of efi'ort has been e xpended 
to un derstand the details of how this occurs (see e.g. lBuddenI 
Il952t) . However, to a good approximation it is enough to set 
the polarisation at the point when tf ~ 2rs to the incoher- 
en t sum of t he polarisation eigenmodes (see the discussion 
in lGinzburjirQTQ.) : 

N = No+Nx 
Nq = ~- cos2x{No - Nx) 
Nu=0 (19) 
Nv = sin 2x(A^o - Nx) 

It is straightforward to show that in terms of the generalised 
Stokes parameters No and Nx are given by (this is true even 
when they are offset by a phase) 



A^o = - (A'' - cos 2xA^Q + sin 2xAv) 
Nx = ^{N + cos 2xNq - sin 2xNv) 



(20) 



Note that, in general, polarisation information will be lost in 
this conversion. This is a refiection of the fact that the space 
spanned by the incoherent sum of the two modes forms a 
subset of the space of unpolarised Stokes parameters. This 
is clear from their respective dimensionalities; the former 
is three dimensional (there are only three degrees of free- 
dom for the decomposition into the two polarisation modes, 
namely their amplitudes and relative phase), while the later 
is four dimensional (/, Q, U, and V, subject only to the 
condition that ^ + + V^). 



4 LOW HARMONIC SYNCHROTRON 

RADIATION INTO COLD PLASMA MODES 

As discussed in the previous section, emission and absorp- 
tion are inherently local processes. As a result it will be 
sufficient in this context to treat them in the LFCR frame, 
and hence in flat space. In this frame it is enough to solve 
the problem in three dimensions and then insert quantities 
in a covariant form. 

Because refractive effects become large only when ui ~ 
ciJ_B,tjp, for there to be significant spectral and polarimetric 
effects it is necessary to have an emission mechanism which 
operates in this frequency regime as well. A plausible can- 
didate is low harmonic synchrotron emission. It is assumed 
that a hot power-law distribution of electrons is responsi- 
ble for the emission while the cold plasma is responsible for 
the remaining plasma effects. In Paper I we did present the 
theory for the warm plasma as well, however, as in the con- 
ventional magnetoionic theory, it is much more cumbersome 
to utilise. 



4.1 Razin Suppression 

A well known plasma effect up on synchrotron emission is 
the Razin suppression (see e.g. iRvbicki fc LightmanllToTSl : 
Bekefi 196(fl . This arises due to the increase in the wave 
phase velocity above the speed of light, preventing electrons 
from maintaining phase with the emitted electromagnetic 
wave, resulting in an exponential suppression of the emission 
below the Razin frequency, 



UJR 



UJp 



LOB 



(21) 



However, as discussed in the Appendix, for the disk model 
we have employed here, typically ujb > i^p and hence the 
Razin effects do not arise. 



4.2 Projection onto Non-Orthogonal Modes 

A significant problem with emission mechanisms in the uj ~ 
lob , Lo p frequency regime is that the modes are no longer 
orthogonal. It is true that for a lossless medium (such as the 
cold plasma), equation |(2J, which defines the polarisation, 
is self-adjoint. However, because of the k'^ differ for the two 
modes, it is a slightly different equation for each mode, and 
hence the polarisations are eigenvectors of slightly different 
hermitian differential operators. In the high frequency limit 
this difference becomes insignificant. 

The energy in the electromagnetic portion of the wave 
(neglecting the plasma portion) is given by 



£ = 



E* 



■e-E _ 1 

4-K An 



47ri 
1 + a 



E 



(22) 



For each mode (Eo and Ex), the dispersion equation gives 



(lo -f A-Kiuicr) ■ Eo,x 



ko,x ® ko,x) ■ Eo,x 



l-k®k -Eo, 



Therefore, with E = E^, 

^=4i2E^|E:.(l-k®kl.E, 



(23) 



(24) 
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However, for a lossless medium it is also true that 



k(g>k -E 



and therefore, 



k]) E* 



1 k( 



. E, = . 



(25) 



(26) 



For a non-degenerate dispersion relation, e.g. that of a mag- 
netoactive plasma, this implies that the the components of 
the polarisation transverse to the direction of propagation 
are orthogonal for the two modes, i.e. 



F* • F,- 



where 



' o.x 



1 - k® k • E, 



E? 



^o,x- (^l-k®kj -Ecx 
As a result it is possible to define £o.x such that 



So,. 



4-K 



and £ = ^ £i . 



(27) 



(28) 



(29) 



i.e. that the electromagnetic energy can be uniquely decom- 
posed into the electromagnetic energy in the two modes. 

Expressions for the Fo,x can be obtained by solving 
for the eigenvectors of the dispersion equation. For the cold 
magnetoactive plasma this gives 



a given polarisation is quoted below: 



o,x 
2 



Fo,x ■ eJ -I- Fo,x ■ ej 



F{x) 



where 



+ [\Fo,x-en\ -|Fo,x-ex| ] G{x) 



2mau 



Sj^eB sin 9 



(33) 



(34) 



/(p) is the distribution function of emitting electron s, Ur is 
the r ay-refractive index (for a suitable definition see iBekefil 
ll96(J) . and F and G have their usual definitions. 



F{x)^x K5{y)dy and G{x) ^ xK2{x) , (35) 

J X 

where the isTs/s and K2/3 are the modified Bessel functions 
of 5/3 and 2/3 order, respectively. The addition factor of 
arises from the difference in the photon phase space, d^k and 
the analogous integral over frequency, Anduj. 

For the adiabatic regime, the emissivities, Jq can 
now be defined: 



Jo,x — 



lpO,X. 

n I ■ 



(36) 



For a power-law distribution of emitting electrons, 
f{p)d?P = C'y~"d'y, this gives 



ko,- 



V2 



l±(l + e)-^/^e|| 



'1/2; 



(30) 



where, (not to be confused with the Levi-Civita pseudo- 
tensor) 



LULUb 



2 cos 9 ujp 



(31) 



9 is the angle between the magnetic field and the wave vec- 
tor, and §11 _x are the fiat space analogues of the basis vectors 
in equation 9 may be defined covariantly by 



cos U — 



{B>^k,) 



B^B, (k^k^+uj^) 



(32) 



Thi s corresponds to the polarisation found in the literature 
(cf.Huddgilliii). 



4.3 Emissivities 



Jo,x 



24tv^ uj2 c{l + s) 

V 4 12 



3 sm I 



' 4 12 



(37) 



The Stokes emissivities and absorption coefficients for an 
emitting hot power law (ignoring effects of order 7"^ as these 
explicitly involve the propagation through the hot electrons) 
are given by 



Jn 
Jq 

lu 



■ Jo + jx 



s 



6 sm ( 



X r 

' 3v = 



+ ^ r £-1 

4 12/ V4 12 



0. 



(38) 

(39) 
(40) 



Note that for low 7 synchrotron can efectively produce cir- 
cular polarisation, namely jy ~ 3/7. The production of cir- 
cular polarisation in this way in environments with large 
Faraday depths will be considered in future publications. 



Because the energy can be uniquely decomposed into the 
energy in each polarisation eigenmode, it is possible to 
calculate the emissivities and absorption coefficients by 
the standard far-field metho d. For synchro tron radiation 
this was originally done by IWestfoldl Il959l The calcula- 
tion is somewhat involved but straightforward and has 
been done in detail in th e subsequent literature (see e.g. 
iRvbicki fc LightmanI 19791) . Consequently, only the result for 
the power emitted (per unit frequency and solid angle) for 



4.4 Absorption Coefficients 

For the adiabatic regime, detailed balance for each mode 
requires that the absorption coefficients are then given by 



oio,x = 



6(jjmc 







6 smy 


V UJ 










1 ± 







'46 



3s + 6 
37+10 



(41) 
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(42) 



In the strongly coupl ed regime, the Stokes absorption 
coefficient matrix (see e.g. iJones fc Q'Deiilll977bl and refer- 
ences therein), 

/ ajv otQ av \ 
aq Qjv 
ajv 
\ Qjv / 

where the Faraday rotation and conversion due to the hot 
electrons have been ignored as a result of the fact that they 
will be negligible in comparison to the Faraday rotation and 
conversion due to the cold electrons. The individual a's can 
be obtained in terms of the ao,x using the fact that the 
energy in the electromagnetic oscillations can be uniquely 
decomposed into contributions from each mode (equation 
Then, 

dN _ dNo_ dNx 
d\ dX dX 

- jo + jx — olqNo — axNx 



(jo +jx) - - {ao +ax)N 



+ ^ cos 2x (qo 



- ax) 



i sin 2x [ao 



ax 



V . 



Therefore, the absorption coefficients may be identified as. 



4.5 



ajv = - [ao + ax) 

aq = -icos2x(Qo - ax) 

av = - sin 2x {ao ~ ax) ■ 



Unpolarised Low Harmonic Synchrotron 
Radiation 



(44) 
(45) 
(46) 



To highlight the role of refraction in the generation of po- 
larisation, an unpolarised emission mechanism is also used. 
To compare with the results of the polarised emission model 
discussed in the previous section, the artificial scenario in 
which the synchrotron emission is split evenly into the two 
modes was chosen. In this case. 



and 



Jo,x — I 



-UP _ - 

Jn — J n ■ 



(47) 



(48) 



with the other Stokes emissivities vanishing. Similarly, the 
absorption coefficients are given by, 

aW^x = aY/' = aN , (49) 
with the other absorption coefficients vanishing as well. 



4.6 Constraints Upon the Emitting Electron 
Fraction 

For refractive plasma effects to impact the spectral and po- 
larimetric properties of an accretion flow, it is necessary that 
it be optically thin. This places a severe constraint upon the 
fraction of hot electrons, / = C/[ne{s — 1)]. In terms of the 
plasma frequency and / the absorptivity is approximately 

' (50) 



aN 



2 

■_ r ^P 

24c'' uj 



CUB . 

3 sm 



With s ~ 2, and ui ~ up, uib, the typical optical depth (not 
to be confused with the affine parameter) is 



hence 



(51) 



where R is the typical disk scale length (here on the order 
of lOM). 



5 RESULTS 
5.1 Disk Model 

Before any quantitative results are presented it is necessary 
to select a specific plasma and magnetic field distribution. 
Here this takes the form of an azimuthally symmetric, thick, 
barotropic disk around a maximally rotating Kerr black hole 
(a ~ 0.98). The magnetic field is chosen to lie upon surfaces 
of constant angular velocity, thus insuring that it does not 
shear. In order to maintain such a field it must also be strong 
enough to suppress the magneto-rotational instability. Fur- 
ther details may be found in the appendix. 



5.2 Ray Trajectories 

Figure shows vertical and horizontal slices of rays propa- 
gated back through the disk discussed in the previous section 
from an observer elevated to 45° above the equatorial plane 
at a frequency cOoc ~ 3ajpmax/4. Note that since the maxi- 
mum occurs at r^q = 2M, the relativistically blue-shifted co 
is approximately 1.8a'pmax placing it comfortably above the 
plasma resonance at all points (assuming Doppler effects do 
not dominate at this point.) 

The refractive effects of the plasma are immediately 
evident with the ex traordinary mode being refra cted more 
so (see discussion in lBroderick fc Blandfordll2003l) . Gravita- 
tional lensing is also shown to be important over a signifi- 
cant range of impact parameters. There will be an azimuthal 
asymmetry in the ray paths due to both the black hole spin 
and the Dopper shift resulting from the rotation of the disk. 
This can be clearly observed in panel (b) Figure 

In panel (a) of Figure the transition between the two 
radiative transfer regimes is also clearly demonstrated. Each 
time a ray passes from the strongly coupled to the adia- 
batic regime it must be reprojected into the two polarisation 
eigenmodes. If the plasma properties (e.g. density, magnetic 
field strength or direction, etc) are not identical to when the 
polarisation had previously frozen out (if at all), this de- 
composition will necessarily be different. As a result, when 
propagating the rays backwards, whenever one passes from 
the adiabatic to the strongly coupled regime, it is necessary 
to follow both polarisation eigenmodes in order to ensure the 
correctness of the radiative transfer. The leads to a doubling 
of the rays at such points. When integrating the radiative 
transfer equations forward along the ray, the net intensity is 
then projected out using equation 12011 . This ray doubling is 
clearly present in panel (a) of Figure where the rays pass 
into the strongly coupled regime and back again as they tra- 
verse the evacuated funnel above and below the black hole. 

Note that the trajectories of the rays depend upon 
iop/ujoo and lub/uoo only (given a specified disk and mag- 
netic field structure, of course), where cuoo is co as measured 
at infinity. Therefore, the paths shown in Figure are valid 
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Figure 1. Shown in panels (a) and (b) are vertical and horizontal cross sections of rays propagating backwards from an observer located 
45° above the equatorial plane. The strongly coupled (adiabatic) regime is denoted by the solid (long-dashed) lines for the ordinary 
(thin) and extraordinary (thick) polarisation eigenmodes. For reference, the null geodesies are drawn in the short dash. In addition, the 
black hole horizon and the boundary of the ergosphere are also shown. 



for any density normalisation of the disk described in the 
appendix as long as ui is adjusted accordingly. 

5.3 Polarisation Maps 

In order to demonstrate the formalism described in this pa- 
per, polarisation maps were computed for the disk model 
described in section IFITl and the appendix^orbiting a max- 
imally rotating black hole as seen by an observer at infinity 
elevated to 45° above the equatorial plane. Each map shows 
Stokes /, Q, U, and V . 

As with the rays trajectories, the particular form of the 
polarisation maps only depend upon a few unitless param- 
eters. These necessarily include uipmax/i^ and uiBmax/^^ as 
these define the ray trajectories. In addition, the relative 
brightness depends upon the optical depth which is pro- 
portional to {ujpjns,x/uj)^{uJBmiix/uj)Mfuj/c. As a result if 
the following dimensionless quantities remain unchanged, 
the polarisation maps shown in the following sections will 
apply (up to a constant scale factor) 

'i'Pmax _ 4 

UJao 3 

WBmax _ 4 

Woo 3 



source on the sky. As a result, Stokes I, Q, U, and V are all 
shown in units of 



M 

T 



2.30 X 10^ 



(52) 



Despite the fact that the form of the polarisation maps 
will remain unchanged if the quantities in equation 1521 re- 
main constant, the normalisation will change by a multi- 
plicative constant in the same way as the source function, 
namely proportional to uj^. However, an additional multi- 
plicative factor arises from the solid angle subtended by the 



M 

— I me LJp max : 



(53) 



where D is the distance to the source. This amounts to plot- 
ting 

^(^)\ (54) 
where Tb is the brightness temperature of the source. 

5.3.1 Unpolarised Emission 

For the purpose of highlighting the role of refractive plasma 
effects in the production of significant quantities of circu- 
lar polarisation. Figure |21 shows Stokes /, Q, U, and V at 
ujoo = 3tJpmax/4, Calculated using the unpolarised emis- 
sion model described in section H31 Immediately noticeable 
are the regions of considerable polarisation surrounding the 
black hole. In addition, the outlines of the evacuated funnel 
above and below the hole are clearly visible. 

Differences in refraction of the two polarisation eigen- 
modes leads two two generic effects: (i) the presence of two 
maxima in the intensity map, each associated with the inten- 
sity maxima in a given polarisation eigenmode; and (ii) a net 
excess of one polarisation, and in particular, circular polari- 
sation. The polarisation changes rapidly at the edges of the 
evacuated funnels because the refraction and mode decom- 
position changes rapidly for modes that just enter the funnel 
and those that pass wide of it. Note that all of the polarisa- 
tion is due entirely to refractive plasma effects in this case. 
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The integrated values for the Stokes parameters are / = 1.3, 
Q = -9.4 X 10"'', [/ = 4.9 X 10-^ and F = 6.2 x 10"^ 
demonstrating that there does indeed exist a significant net 
circular polarisation. 

Figure |21 may be compared with Figure |3] in which 
Stokes I, Q, U, and V are shown at ujca = Slop max for the 
same unpolarised emission model. In the latter case the re- 
fractive effects are significantly repressed. This demonstrates 
the particularly limited nature of the frequency regime in 
which these types effects can be expected to occur. In 
this case there still does exist a net circular polarisation, 
now with integrated values I = 1.0, Q — —4.8 x 10"'', 
[/ = 2.4 X 10"^ and V = 1.2 x 10"^ 

5.3.2 Polarised Emission 

In general, synchrotron emission will be polarised. As a re- 
sult it is necessary to produce polarisation maps using the 
emission model described in sections 14.31 and 14.41 In this 
case a net polarisation will exist even in the absence of any 
refraction. In order to compare the amount of polarisation 
generated by refractive effects to that created intrinsically. 
Figure 2] shows Stokes /, Q, U, and V calculated using the 
polarised emission model and ignoring refraction (i.e. set- 
ting the rays to be null geodesies) for uJoo = 3a'p,max/4. 
Strictly speaking, this is a substantial over estimate of the 
polarisation. This is because, in the absence of refraction, in 
principle it is necessary to include Faraday rotation and con- 
version in the transfer effects considered. As a result of the 
high plasma density and magnetic field strengths, the Fara- 
day rotation and conversion depths for this system should be 
tremendous for non-refractive rays, effectively depolarising 
any emission. 

In comparison to Figures |21 and |21 the general morphol- 
ogy of the polarisation maps are substantially different. In 
addition, the amount of linear polarisation is significantly 
larger, having an integrated value of over 60% compared 
to less than 0.1% in Figure |5| and less than 10~^% in Fig- 
ure El This calculat ion can be compared to that done by 
iBromlev et"ai]|200ll . In both it was assumed that the rays 
were null geodesies . In both Faraday r otation /conversion 
were neglected (in iBromlev et alJ l200ll because for their 
disk model it was a ssumed to be negligible.) However, in 
IBromlev et"ai]l200ll it was also assumed that the radiative 
transfer could always be done in the adiabatic regime. As a 
result, the net polarisation was determined entirely by the 
emission mechanism. However, as discussed in section 13.11 
this is only possible in the strongly coupled regime. In this 
case, the dichroic terms in equation (1421 provide the source 
of circular polarisation, even in the absence of a circularly 
polarised emission, resulting from the different absorption 
properties of the two polarisation eigenmodes. This is what 
leads t o the presence of circ ular polarisation in Figure 2] but 
not in IBromlev et"alll200ll . In this case, the integrated val- 
ues of the Stokes parameters are / = 1.1, Q = 6.0 x 10^^, 
U = -4.9 X 10"^, and V = 6.9 x . The vertical feature 
directly above the black hole in panels (b) and (c) are asso- 
ciated with the rapid decrease in the magnetic field strength 
in the evacuted funnel above and below the black hole and 
are due to the geometric transfer effect discussed in section 

lO 

Finally, in Figure |S| both refractive effects and the po- 



larised emission mechanism are included (again at ujoo = 
3ti;p,niax/4). Many of the qualitative features of Figure (21 
still persist. The integrated values of the Stokes parame- 
ters are I = 1.3, Q = -2.2 x 10"^ U = 1.2 x 10"", and 
V = 1.4 X 10~^. While the intrinsic polarisation in the emis- 
sion does make a quantitative difference, it is clear that in 
this case the generic polarimetric properties are dominated 
by the refractive properties. This is most clearly demon- 
strated by noting the strong supression of linear polarisa- 
tion. In Figure |K|the linear polarisation fraction is less than 
0.2% as compared with nearly 60% in Figure 2] 

5.4 Integrated Polarisations 

Figure |S| shows the Stokes parameters as a function of fre- 
quency for when only polarised emission is considered, only 
refractive plasma effects are considered, and when both are 
considered. There are two notable effects due to refraction: 

(i) the significant suppression of the linear polarisation, and 

(ii) the large amplification of circular polarisation. The lin- 
ear polarisation is decreased by at least two orders of magni- 
tude, and in particular, at least two orders of magnitude less 
than the final circular polarisation. On the other hand, the 
circular polarisation is more than doubled at its peak, and 
increases by many orders of magnitude at higher frequencies. 
Nonetheless, by cOoo = lOtjpmax, both polarisations are less 
than one tenth of their maxima. As a result, it is clear that 
this mechanism is restricted to approximately one decade in 
frequency, centred about ojpmax. 

Figure [7| shows the circular polarisation fraction as a 
function of frequency for the same set of cases that were 
depicted in the previous figure. As can be seen in Figure [S] 
the circular and linear polarisation spectral index are ap- 
proximately equal, and both are softer than that of the to- 
tal intensity. The result is a decreasing circular polarisation 
fraction with increasing frequency. 



6 CONCLUSIONS 

We have presented refraction as a mechanism for the gen- 
eration of polarisation when cuoo ~ wp, wp. That this will 
typically result in mostly circular polarisation is a result of 
the fact that the polarisation eigenmodes are significantly 
elliptical only when the wave- vector and the magnetic field 
are within oip/o) of perpendicular, which is usually a small 
number near the surface where the polarisation freezes out. 
In addition to producing circular polarisation, this mech- 
anism also significantly suppresses linear polarisation. Be- 
cause it does require significant refraction to take place, it is 
necessarily limited to approximately a decade in frequency, 
making it simple to identify. 

As shown in section 1^^ the resulting circular polarisa- 
tion will be softer than the intensity. However, because of op- 
tical depth effects, as the observation frequency increases the 
polarimetric properties of such a system will be dominated 
be increasingly smaller areas. As a result, the fractional vari- 
ability in such a system would be expected to increase with 
frequency. Furthermore, even though the emission may arise 
from a large region, the polarimetric properties will continue 
to be determined by this compact area, making it possible 
to have variability on time scales short in comparison to 
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Figure 2. Stokes I, Q, U, and V per unit A/^ are shown in panels (a), (b), (c), and (d), respectively, for the unpolarised emission 
mechanism described in section 14.51 and the disk model described in section 15.11 and appendix ^ orbiting a maximally rotating black 
hole from a vantage point 45° above the equatorial plane at the frequency uicxi = 3ajp inax/4. The contour levels are at 0.2 (dashed) 
and 0.6 (solid) of the maximum values shown on the associated colorbars. The integrated fluxes over the region shown are / = 1.3, 
Q = -9.4 X 10"'', U = 4.9 X 10~^, and V = 6.2 x 10~^. All fluxes are in units of (M/ D)~^uij, as discussed above equation l53t . 



those associated with the emission region. In addition, vari- 
ability in the circular polarisation would be expected to be 
correlated with variability in the integrated intensity at fre- 
quencies where the emission is dominated by contributions 
from close to the horizon (e.g. X-rays). 

Possible applications to known astrophysical sources in- 
clude the Galactic Centre (at submm wavelengths) and ex- 
tinct high mass X-ray binaries (in the infrared). These will 
be disscussed in further detail in an upcoming paper. 
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APPENDIX A: A THICK DISK MODEL 

In general, the innermost portions of the accretion flow will 
take the form of a thick disk. The equation for hydrostatic 
equilibrium in the limit that is given by 



= -d^ lnE + 



1 - flL 



(Al) 



where here F is t he adiabatic index, E = —Ut. ~ u"^ /u^ , 
and L — —u^/ut iBlandford fc Be gelman"2003^ . Note that, 
given the metric, any two of the quantities E, Q., or L, may 
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Figure 3. Same as Figure 1^ except with u)oo = Scjpmax- The integrated fluxes over the region shown are I = 1.0, Q = —4.8 X 10 ^, 
1/ = 2.4 X 10^'^, and V = 1.2 X lO^''. All fluxes are in units ol {^^I/ D)^^ujp as discussed above equation I5:H . 



be derived from the third. Explicitly, Q and L are related 
by 



gtt _|_ gt,Pl 



(A2) 



and the condition that u'^u^ — u^ut + u'^u^ = — 1 gives E 
in terms of fl and L to be 



E = 



[- [g" + g'^L){l-nL)] 



-1/2 



(A3) 



In principle this should be combined with a torque bal- 
ance equation which explicitly includes the mechanism for 
angular momentum transport through the disk. However, 
given a relationship between any two of the quantities E, Q, 
and L specifies this automatically. Thus the problem can be 
significantly simplified if such a relationship can be obtained, 
presumably from the current MHD disk simulations. 



Al Barotropic Disks 

For a barotropic disk the left side of equation I|A1^ can be 
explicitly integrated to define a function H: 



H ■ 



dP 



(A4) 



which may be explicitly integrated for gases with constant 
r to yield 



H = \n[ 1 



r P 

Therefore, reorganising equation lAl^ gives 

Q.du.L 



d^{H-lnE) 



1 - fii ' 



(A5) 



(A6) 



which in turn implies that f2 is a function of L alone. Speci- 
fying this function allows the definition of another function 
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Figure 4. Same as Figure [^except using the polarised emission mechanism (described in sections 14. 31 and 14. 41 and ignoring refractive 
plasma effects. The integrated fluxes over the region shown are / = 1.1, Q = 6.0 X 10~^, U = —4.9 X 10~^, and V = 6.9 X 10~^. All 
fluxes are in units of {!^/ D)~'^ujp as discussed above equation I53i . 



n{L)dL 
1 - n{L)L 



(A7) 



Using their definitions, it is possible to solve Q — n{L) for 
L{x'^) and hence H(a;''). Then H and H are related by. 



H = Ho + lnE -E, 



(A8) 



which then may be inverted to yield p{Ho — In i? + H) . In- 
verting H for p then yields p{x'^). The quantity Hq sets the 
density scale and may itself be set by choosing p at some 
point: 



Ho = H{po)-iliiE^S)ix^). 



(A9) 



A 1.1 Keplerian Disk 

As a simple, but artificial, example of the procedure, a Ke- 
plerian disk is briefly considered in the limit of weak gravi- 
tating Schwarzschild black hole (i.e. r 3> M). Note that this 
cannot be done in flat space because in equation ijAlfl the 
gravitational terms are present in the curvature only. For a 
Keplerian flow, = ^M/{rsme)^ ~ M^L~^. In that case 
using the deflnition of H gives 



dL 



d£ 



(AlO) 



where £ = L/M. However, In E is given by 



\nE - 



In = In \/l - ^ - In vT^^^^, 

(All) 
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Figure 5. Same as Figure |4] except including refractive plasma effects. The integrated fluxes over the region shown are / = 1.3, 



Q = -2.2 X 10-3, C/ = 1.2 X 10"*, and V = 1.4 X 10"^. All fluxes are in units of (M/D)- 



as discussed above equation 1531 . 



and hence, 



H = Ho~lnE + : 



2M 

7?o-ln\/l- — +ln(l 



Ho + 



M 



M 
r sin 9 



(A12) 



where i — ^ r sm6/M and the weakly gravitating condition 
were used. As expected, along the equatorial plane H , and 
therefore p, is constant. For points outside of the equato- 
rial plane pressure gradients are required to maintain hy- 
drostatic balance. 



A 1.2 Pressure Supported Disk 

Accretion disks will in general have radial as well as vertical 
pressure gradients. Inward pressure gradients can support 
a stable disk inbetween the innermost stable orbit and the 



photon orbits, thus decreasing the radius of the inner edge of 
the disk. Around a Schwarzschild black hole this can bring 
the inner edge of the disk down to 3Af. In a maximally 
rotating Kerr spacetime this can allow the disk to extend 
down nearly to the horizon. 

Far from the hole, accreting matter will create outward 
pressure gradients. An angular momentum profile appropri- 
ate for a Kerr hole which goes from being super to sub- 
Keplerian is 



i(rcq) = < 



if Tcq < Til, 



otherwise 



gtt _|_ gt4,l 



(A13) 
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Figure 6. The log of the integrated intensity, total linear polar- 
isation, and circular polarisation are shown as a function of the 
observation frequency at infinity for when only polarised emission 
is considered (open triangles), only refractive plasma effects are 
considered (open squares), and when both are considered (filled 
circles). As in Figures the disk model described in section 
15. H and appendix ^ orbiting a maximally rotating black hole is 
viewed from a vantage point 45° above the equatorial plane. All 



fluxes are in units of (M/D) 
tion 1531 . 



as discussed above equa- 



where both L and Q, are parametrised in terms of the equato- 
rial radius, rcq. The condition that L reduces to the angular 
momentum profile of a Keplerian disk for radii less than 
the inner radius ensures that no pathological disk struc- 
tures are created within the photon orbit. The constants 
ci and C2 are defined by the requirement that at the inner 
edge of the disk, rinnor, and at the density maximum, rmax, 
the angular momentum must equal that of the Keplerian 
disk. In contrast, lo is chosen to fix the large r behaviour 
of the disk. The values chosen here were rinner ~ 1.3M, 
''max = 2M, and lo = 0.1. The value of Hq was set so 



that H{r 



lOOM) = 0, thus making the disk extend 



In addition to defining Q and L it is necessary to define 
P{p)- Because the gas in this portion of the accretion fiow is 
expected to be unable to efficiently cool, F = 5/3 was chosen. 
The proportionality constant in the polytropic equation of 
state, K, is set by enforcing the ideal gas law for a given 
temperature (To) at a given density (po)- Thus, 



P{p) = po 



kTo 



5/3 



(A14) 



Note that po and To provide a density and temperature scale. 
A disk solution obtained for a given po and To may be used to 
generate a disk solution for a different set of scales simply 
by multiplying the density everywhere by the appropriate 
constant factor. 



Figure 7. Shown is the circular polarisation fraction as a function 
of the observation frequency at infinity for when only polarised 
emission is considered (open triangles), only refractive plasma ef- 
fects are considered (open squares), and when both are considered 
(filled circles). As in Figures the disk model described in sec- 
tion lS.lj and appendix^orbiting a maximally rotating black hole 
is viewed from a vantage point 45° above the equatorial plane. 



A2 Non-Sheared Magnetic Field Geometries 

The disk model discussed thus far is purely hydrodynamic. 
Typically, magnetic fields will also be present. In general, it 
is necessary to perform a full MHD calculation in order to 
self-consistently determine both the plasma and magnetic 
field structure. However, an approximate steady state mag- 
netic field can be constructed by requiring that the field lines 
are not sheared. 

To investigate the shearing between two nearby, space- 
like separated points in the plasma, x'^ and xi^, consider the 
invariant interval between them: 



As = Ax'^Axn where Ax^ = xii 



(A15) 



The condition that this doesn't change in the LFCR frame 
is equivalent to 



dAs^ 
ds 



0. 



(A16) 



Expanding in terms of the definition of As gives. 



d dx" 
—g^r.Ax^Ax" = g^^.a — Ax^Ax" 
ds ' ds 



+ 2g^^Ax' 



, dAx" 
ds 



0. (A17) 



Note that by definition. 



dx^" _„ , dAx" 
—— = u'^ and — ; — 

ds ds 



u „Ax 



(A18) 
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Figure Al. Shown are the contours of the density and azimuthal 
velocity as measured by the zero angular momentum observer, 
and the magnetic field lines. Starting at the density maximum 
{''eq = 2M and z = 0), the density is contoured at levels 10""'^ 
to 10-*-5 times the maximum density in multiples of 10 ^ . From 
left to right, the velocity is contoured at levels 2~^'^c to 2~^c 
in multiples of 2~'^'^. In order to provide a distinction between 
the velocity contours and the magnetic field lines, the velocity 
contours are terminated at the disks surface. 



Hence, 



ds 



= {gfj.u.irlf + 2u^,i, — 2g^a,iAf) Ax'^ Ax" 

= 2 - Vl^Ucr) Ax^'Ax" 

= 2 (V„u,.) Ax^'Ax" = . (A19) 



The final equality is easy to understand from a geomet- 
rical viewpoint; for there to be no shearing, there can be 
no change in the direction of Aa;*^ of the component of the 
plasma four- velocity along Ax'^ . 

That a steady state, axially symmetric magnetic field 
must lie upon the non-shearing surfaces can be seen di- 
rectly by considering the covariant form of Maxwell's equa- 
tions. In particular V^^'"' = 0, where *F'^'^ is the dual of 
the electromagnetic field tensor, which in the absence of 
an electric field in the frame of the plasma takes the form 
= B^'Tt - B"u''. Therefore, 

- B^u^V^B" - B^.B^V.u" 
= -B^B^V.u^ = , (A20) 

where the first three terms vanish due to the symmetries 
and the requirement that B'^u^ = 0. This is precisely the 
non-shearing condition obtained in equation (IA19II . 

For plasma flows that are directed along the Killing 



where t'^ is the time-like Killing vector, it is possible to sim- 
plify the no-shear condition considerably. 

Ax^Ax^-^^jf 

= Ax^Ax^ (^u'Vf.f + J2 « V^C>" 

-I- Ax'^Ax^ ^"dfiU + ^2 &df^u 
= AxtAx^d^u + Ax.Ax^d^u = , (A22) 

i 

Where terms in the first parentheses vanish due to Killing's 
equation. The additional constraint that Ax^u'^ — gives 



Axt = - ^ i^iAxi , 



(A23) 



where Qi = u ju is a generalisation of the definition of S7 
at the beginning of the section. Inserting this into equation 
<A22l l and simplifying yields 



^ Ax.Ax'^d^Q., 



0, 



(A24) 



i.e. the no shear hypersurfaces are those upon which all of 
the Sli are constant. 

For the plasma fiows considered in ilAll the plasma ve- 
locity is in the form of equation IIA2H where the space-like 
Killing vector is that associated with the axial symmetry, 
0''. Thus with f20 = fi, the no-shear condition for this class 
of plasma flows is 



Ax^'dM = . 



(A25) 



Note that while we have been considering only axially sym- 
metric plasma flows, this no shear condition is more gener- 
ally valid, extending to the case where f2 is a function of t 
and as well as r and Q. However, in this case it is not the 
perfect-MHD limit of Maxwell's equations. 

For a cylindrically symmetric disk, the no-shear condi- 
tion may be used to explicitly construct the non-shearing 
poloidal magnetic fields by setting 



e'' = BQ.,e and B" 



(A26) 



Once the magnitude of B^ is determined at some point along 
each non-shearing surfaces (e.g. in the equatorial plane), it 
may be set everywhere by V fj.B'^ — B^TtV vUf_i — 0, which 
comes directly from Maxwell's equations in covariant form 
and B'^u^ — 0. Inserting the form in equation <A26t into 
the first term gives 

V9 



1 

1 

71 

B^d^ In ^B 



(A27) 
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The second term can be simplified using equation (IA21L 

+ B^''u''"V f^Uv 
= B^ {utd^u +u^d^uQ.) 
— C (Tit + itu^) df_tU + u^uB'^df_tVl 



-B^'dalnt 



(A28) 



where the stationarity and axially symmetry have been used 
in the third step and the no-shear condition was used in the 
final step. Therefore, the magnitude B can be determined 
by 



V^S" - B'^lfV^Uf, = B^d^ In ^B - B^d^ Intt* 

= B^d^ In ^ = , (A29) 
u 

and hence 



constant 



(A30) 



along the non-shearing surfaces. If B is given along a curve 
which passes through all of the non-shearing surfaces (e.g. in 
the equatorial plane), B'^ is defined everywhere through 
equations IIA26t and 1A30II . 



A2.1 Non-Shearing Magnetic Fields in a Cylindrical Flow 

An example application of this formalism is a cylindrical 
flow in fiat space. In this case, f2 is a function of the cylin- 
drical radius zu = rsin^. The Keplerian disk is a specific 
example with D, = tz!~^^'^. The direction of the magnetic 
field is determined by, 



il,r = — sm t 
dzu 



and Q^e ~ -f-r cos9 . (A31) 
aw 



The magnitude, B is given by 



Vl-r-^sin^ 60.'^ 



and thus 



B = -b{zu) 



(A32) 



(A33) 



where the particular form of 6(1*7) depends upon the partic- 
ular form of /(fi). Therefore, 



B^ — h{vj) cos 9 and 



B" 



-h{vj)-sme, (A34) 



which is precisely the form of a cylindrically symmetric ver- 
tical magnetic field. 



A2.2 Stability to the Magneto -Rotational Instability 

A sufficiently strong non-shearing magnetic field configura- 
tion will remain stable to the magneto-rotational instability 
(MRI). The criterion for instability to the MRI is 



(k • va) < -r 



dr 



(A35) 



where k is the wave vector of the unstable mode and va 
is the Alfven velocity ( Hawlov & Balbus'lOOS"). For a nearly 
vertical magnetic field geometry, stability will be maintained 
if modes with wavelength less than twice the disk height (h) 
are not unstable. With 



B 



VA 



a Keplerian disk will be stable if 



(A36) 



(A37) 



A conservative criterion may be obtained by approximating 
h ~ hor for some constant of proportionality ho, hence 



LUB 



M 
r 



0.3, 



(A38) 



for ho ~ 0.1 and r ~ 7 which are typical for the disk pictured 
m Figure IXn 

Comparison to equipartion fields can provide some in- 
sight into how unrestrictive the stability criterion really is. 
Given /3 = Pgas/^mag and the ideal gas law it is straight 
forward to show that 



uip 



2kT 



(A39) 



where T is the ion temperature. Because the ion temperature 
in a thick disk will typically be on the order of or exceed 10^^ 
K, the equipartition ujb {P — 1) will be at least an order of 
magnitude larger than ujp. As a result the field needed to 
stabilise the disk against the MRI is an order of magnitude 
less than equipartition strength, and hence is not physically 
unreasonable. 



A2.3 Magnetic Field Model 

Considering the restriction placed upon the magnetic field 
strength discussed in the previous sections, B was set such 
that in the equatorial plane 

u}B = iop + ri{r + 10M)~^^* , (A40) 

where the second term provides a canonical scaling at large 
radii. Here rj was chosen to be 0.01. 
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